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Abstract. Within a test-particle description of the acceleration process in parallel nonrelativistic 
shocks, we present an analytic treatment of the electron injection. We estimate the velocity distribution 
of the injected electrons as the product of the post-shock thermal distribution of electrons times the 
probability for electrons with a given velocity to be accelerated; the injection efficiency is then evaluated 
as the integral of this velocity distribution. We estimate the probability of a particle to be injected as that 
of going back to the upstream region at least once. This is the product of the probability of returning 
to the shock from downstream times that of recrossing the shock from downstream to upstream. The 
latter probability is expected to be sensitive to details of the process of electron thermalisation within the 
(coUisionless) shock, a process that is poorly known. In order to include this effect, for our treatment we use 
results of a numerical, fully kinetic study, by Bykov & Uvarov (1999). According to them, the probability of 
recrossing depends on physics of thermalisation through a single free parameter (F), which can be expressed 
as a function of the Mach number of the shock, of the level of electron-ion equilibration, as well as of the 
spectrum of turbulence. It becomes apparent, from our analysis, that the injection efficiency is related to 
the post-shock electron temperature, and that it results from the balance between two competing effects: 
the higher the electron temperature, the higher the fraction of downstream electrons with enough velocity 
to return to the shock and thus to be ready to cross the shock from downstream to upstream; at the same 
time, however, the higher the turbulence, which would hinder the crossing. 

Keywords: Shock waves, acceleration of particles, cosmic rays, supernova remnants 

PACS numbers: 98.38.Mz, 95.30.Lz 

1 Introduction 

Strong coUisionless shocks are present in various astrophysical objects, and under a wide range 
of conditions. These shocks effectively heat the gas and are also beheved to accelerate a fraction 
of particles up to very high energies. Analyses based on multi-frequency observations allow one 
to determine properties of both thermal and nonthermal component. In the case of supernova 
remnants (hereafter SNRs), optical and UV emission are typically thermal, radio is nonthermal, 
while thermal and nonthermal emission may coexist in X-rays. The momentum distribution of the 
accelerated particles is locally well approximated by a power law; this can be inferred from the 
power-law synchrotron spectra, in the case of electrons; while in the case of ions it can be measured 
directly in the energy distribution of cosmic rays (if one accepts the "SNR Paradigm" for the origin 
of galactic cosmic rays). 

Diffusive acceleration (Fermi acceleration) is beheved to be the dominant process that allows 
particles to gain energies in excess of typically thermal values. The standard theory of diffusive 
acceleration, in test-particle approximation (see e.g. the review of Jones & EUison 221), shows 
that a power-law distribution develops at high energies. The spectral index of this high-energy 
population depends on the shock compression ratio, while it does not depend at all on the original 
energy distribution of the injected particles. In fact, a simplified way to describe the overall particle 
evolution, from nearly thermal to very high velocities, is to treat particle injection and acceleration 
as two separate problems. The injection problem consists in finding out the initial momentum 
distribution of that fraction of (originally thermal) particles that can enter the acceleration process, 
i.e. to make at least one acceleration cycle. The acceleration problem, instead, consists in following 
the evolution of the distribution of these particles along all next acceleration cycles. 

In order to keep the number of free parameters low, while modelling more effectively the emis- 
sion in all observed spectral ranges, one needs to introduce a physically self-consistent scenario 



for the thermal and nonthermal populations. For instance, the standard model of particle accel- 
eration constrains the slope of the electron distribution at high velocities, but does not predict 
its normalization: in other terms, the injection efficiency (i.e. the fraction of particles that enter 
the acceleration process) is poorly known, because this process is sensitive to physical details not 
included in the standard model of Fermi acceleration. The level of electron-ion equilibration or, 
alternatively, the electron temperature is another key quantity hard to determine "a priori" in col- 
lisionless shocks. While in models of SNR shocks injection and equilibration efficiencies are taken 
as independent free parameters, in the reality both depend on the physical conditions within the 
shock transition, and therefore they are not independent. Goal of this paper is in fact to investi- 
gate, in nonrelativistic SNR shocks, a possible connection between electron injection and thermal 
equilibration. 

A self-consistent treatment of injection and acceleration must include a microphysical model of 
particle-wave interactions in the plasma. A few physical processes have been proposed to account 
for the electron injection (see Malkov & Drury p^j for a review). The scattering of electrons is 
suggested to be due to some ion-generated instabilities (Bykov & Uvarov 0; hereafter BU99), 
whistler waves QJ] and lower-hybrid waves from ions j^I- These models mainly address the 
plasma microphysics; while only BU99, to our knowledge, are able to model the formation of the 
post-shock electron distribution. 

In this paper, we approach the problem in a simplified way. We assume the presence of 
scattering centres, without concentrating on their nature, but simply assuming that they match the 
following requirements: i) the interaction with these scattering centres generates a nearly isotropic, 
Maxwellian velocity distribution of particles on timescales not longer than one collision time; ii) 
the timescales for (wave-mediated) isotropization and energy exchange between electrons are both 
smaller than the (wave-mediated) electron-ion equihbration time; iii) the scattering centres play 
at the same time the role of thermaHsing, within the shock, the incoming particle population to 
the post-shock temperature and that of driving the process of diffusive acceleration. When we will 
need to use more specific properties of the wave-particle interaction (in Sect. 3.2), we will refer to 
the results of BU99 on the electron kinetics. 

There are in general three ways to calculate the post-shock momentum distribution of particles: 
either by solving the kinetic equations, or by making a hybrid simulation (which is however unable 
to model the momentum distribution of electrons because electrons are treated as a fluid), or finally 
by extending the individual particle approach of Bell jl]. He has estimated the probability for a 
particle to return to the shock from downstream and has shown that in this way one obtains a 
power-law distribution for the accelerated particles with velocities v ^ Vs (where K, is the shock 
velocity). 

In the present paper we propose to extend the Bell approach to the problem of injection by 
introducing the probability to recross the shock from downstream to upstream. This probability 
is connected to the process of thermalisation of the incoming flow within the shock. This fact has 
been shown by Malkov ^Hl in the case of protons. Namely, the idea that ions are prevented from 
backstreaming by the self-generated waves (which also participate in thermaHsation of ions) has 
allowed Malkov to obtain an analytic solution of the injection problem for protons. The main point 
in his thermal leakage theory is that only those protons that can "leak" upstream are injected into 
the Fermi process. 

We use the same idea in our approach to electrons, even though we are not tight to any 
specific kind of interaction. We consider only the case of parallel shocks, namely when the ambient 
magnetic field is parallel to the shock normal. 

The plan of the present paper is as follows. In Sect. El we deal with the injection problem by 
developing an analytic treatment for calculating the injection efficiency as well as for determining 
of the distribution of particles which are able to be accelerated. Sect, ideals with the process of 
thermalisation and its influence on injection and, using results from BU99 model, give quantitative 
estimations on the injection efficiency. Sect. concludes. 



2 Injection efficiency and initial distribution 



2.1 Efficiency of electron injection 

Let us assume that all electrons are injected into the acceleration process from the downstream 
thermal population, i.e. we do not invoke seed particles with velocities already much higher than 
the thermal velocity. Their distribution is then well approximated by nesfu, where 

/m(?/) = 4=2/^exp(-j/2) (1) 

is a normalized Maxwellian, isotropic in the fluid comoving frame. We have introduced the reduced 
momentum y = p/pth, which is also equal to the reduced velocity v/vth, as long as non-relativistic 
particles are considered, as it is the case at injection. Thermal momentum and velocity arc defined 
by Pth — 'TioWth ~ ^/^rrickTcs, where Tcs is the post-shock electron temperature. We consider a 
fully ionized H+He gas (with nne = O.lnn, for a moan mass per particle //, = 0.609), and a strong, 
unmodified shock. For an adiabatic index 7 = 5/3, the shock compression ratio is cr = 4 but, for 
the sake of generality, in the following formulae we shall allow for a general a. Therefore, the ratio 
of the electron thermal velocity, Vth, to Vg is 



^= /^^, (2) 

where Xo = We/(/xmp) ~ 8.94 x 10""*. The factor Xs = Tes/Tg, where Tg is the mean shock 
temperature, accounts for the thermal equilibration level between electrons and ions immediately 
after the shock, and ranges from Xo (no equilibration) to 1 (full equilibration). 

Introducing the simple-minded assumption that only particles in the high-velocity tail of the 
MaxweUian distribution can be accelerated, it is easy to Hnk the minimum momentum of this tail, 
Ptaii = ytaiiPth, to the injection efficiency <; (i.e. the fraction of accelerated particles). One has 
just to solve the equation fuiv) = which gives for instance ytaii = 2.85 for <r = 10"'^ and 
Vtaii = 3.91 for <; = 10~^. It is worth noticing that, for reasonable values of <;, this integral is 
dominated by particles with y ~ ytaii, with j/taii of order of unity: it is apparent from this example 
that injection involves mostly particles with velocities of the order of the thermal one, and not 
only those with v ^ Vg. 

In the above estimation, we have assumed that all particles with y > yma, and only them, are 
accelerated. In order to find out the injection efficiency in a more general case, we introduce the 
probability V{y) for a particle with velocity v = yvth to be accelerated, i.e. to recross the shock 
from downstream to upstream at least once. This probability yields the fraction of particles, with 
a given velocity, which can be accelerated; while the Maxwellian distribution in turn gives the 
number density of particles with that velocity. Thus, for an isotropic velocity distribution, the 
fraction of accelerated particles (injection efficiency) is given by the integral 

00 

? = J ■P{y)fM{y)dy. (3) 


In other terms, the distribution of particles injected into the acceleration process is 

fini{y) = V{y)fM{y). (4) 

The probability P(?y) in turn can be estimated as the product of the probability, P^, that a particle 
returns to the shock from downstream, times the probability, Pc, that this particle crosses the shock 
moving upstream. The next two subsections will be devoted to estimate these two probabilities. 

We wish to point out that a common misconception lies underneath the Fermi acceleration 
approach, namely that the electrons must enter this process having already a velocity much higher 
than Vs. This is usually obtained, by requiring either i) that the electron temperature is close 
to equipartition, or ii) that only electrons in the high-energy tail of the Maxwellian distribution 
enter into the acceleration process, or finally iii) that some unknown pre-acceleration mechanism 
takes place to accelerate electrons to the required velocity regime. The condition v ^ 14 is in fact 



very useful to simplify the mathematical treatment of the process, but in our belief is not strictly 
required by physical arguments. In the present paper we will show instead i) that electrons may 
be injected efficiently also when their temperatures is far from equipartition, ii) that, in order 
to have reasonably high injection efficiencies, Xs has to be considerably less than unity; in other 
words, the velocities of the majority of the injected electrons must not be not too far from the 
thermal velocity and the minimum injection momentum can even be much smaller than thermal 
one, and finally iii) that there is no physical need for an independent pre-acceleration process, if 
the treatment of the acceleration is modified in order to account also for relatively low particle 
velocities (this can be done by introducing the probability of crossing the shock. Sect. 12.311 . 



2.2 Probability of returning to the shock 



In the case of isotropic velocity distribution in the downstream flow, the probability for particles 
with velocity v to return to the shock from downstream is given by the ratio of the upstream and 
downstream fluxes |12) : 



J {U2 + Vx) dVx 
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(5) 



where: U2 = Ks/cr is the velocity of the downstream flow, in the shock reference frame; v is the 
velocity, in the downstream flow reference frame, of the particle that has just reached the shock; 
H(w — U2) is the Heaviside step function (meaning that U2 is the minimum value of v that allows a 
particle to return to the shock). In general, in the paper we label quantities refering to upstream 
with "1" and quantities refering to downstream with "2". Coordinates are defined in such a way 
that, in the reference frame of the shock, the flow moves along the a;-axis in the positive direction. 

It is useful to re-write Eq. Q to fix a lower boundary to the reduced momenta of the injected 
particles 

1/2 
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(6) 



The quantity j/min is always less than (2(ct — 1))^^/'^ (i.e. Pmin < 0.4pth for a — A) and can be 
much smaller than that if Xs ~ 1 (which means Vt\i ^Vs). This means that, the higher the level of 
electron-ion equilibration, the higher the electron thermal velocity compared to V^, and thus the 
higher the fraction of electrons able to return to the shock from downstream. 



2.3 Probability of crossing the shock 

The standard theory of diffusive acceleration ^ impHcitely assumes Pc — 1 , which means that the 
particle mean free path A is longer than the thickness. Ax, of the shock transition region. This 
condition applies only for particles with high enough velocity (v ^ V"s). 

On the contrary, the evolution of particles with lower velocities is affected by scatterings 
within the shock transition. In fact, the mere existence of a shock implies that the incoming am- 
bient plasma must be thermaHsed, within the shock transition region, by some kind of scatterings 
centres. Also particles that enter the shock transition region from downstream, as long as they 
have velocities similar to thermal particles, must experience a similar rate of scatterings. Thus, 
also for them A < Ax. 

In the presence of scatterings, only a fraction of these particles will succeed crossing the shock 
and finally reaching the upstream region. In general, modelling this process is very complex. Here 
we will present a simplified treatment, based on some approximations. The first of them is diffusive 
approximation, which requires that mean free paths are smaller than the shock thickness, and that 
the velocity distribution is nearly isotropic. 

However, this assumption is invalid near the downstream boundary of the shock layer. In fact, 
the original distribution of downstream particles which return to the shock is highly anisotropic, 
since all particles entering the shock have, in the shock reference frame, an x-component opposite 
to the flow velocity. The above assumption is anyway valid over most of the volume, provided 



that isotropization processes within the shock are very efficient. Namely, we require that the 
length scale for isotropization is of the order of one mean free path (similarly to what happens for 
Coulomb collisions between similar particles). 

The estimation of Pc{v) is generally very complex. Here we use a crude approximation (based 
on the so-called "modulation" equation, see e.g. |^) and write 

Pc = exp{-{u)Ax/K). (7) 

where (u) « {ui + U2)/2 = U2(cr + l)/2 and k is the diffusion coefficient. 

We want to point out that probability of crossing is closely related to the thermalization level 
Xs- The thickness Aa; may be derived from the condition that the temperature of the incoming fluid 
increases to the post-shock value Tes = Xs^s while the fluid moves through the shock transition 
from upstream to downstream. In this way, the two problems - injection and thermalisation - 
become closely connected. This can be seen by rewriting 

Ax = {u)Ah2 = {u) I (^^ j dx (8) 



where Ati2 is the time it takes to a fluid element to cross the shock, moving from its upstream 
boundary to the downstream one. In general, it is necessary to introduce a microphysical model of 
thermalisation in order to obtain explicitely the functional dependence of the rate of thermalization 
dx/dt. 

For the diffusion coefficient we use the standard formula k — A'u'/3, where v' is the velocity of 
a particle in the local reference frame of the flow and A' is the particle mean free path with respect 
to scatterings within the shock transition. A further assumption behind this formula is that the 
scattering centres are frozen into the fluid. This is, for instance, the case in the BU99 model for 
the electron kinetics in a strong shock. In this model, particles are scattered by the ion-generated 
Alfvenic waves, and the Alfvenic speed is much lower than the shock velocity. In case of electron 
diffusion in presence of magnetic field turbulence, it is common to parametrize the mean free path 
as A' = rjrg, where rg = p'c/eB is the gyroradius and rj accounts for the level of turbulence. We 
concentrate here on particles with velocities not much larger than Vs , and therefore we will use the 
nonrelativistic formula for A'. For such parameterization of the mean free path. A' may be written 
as A' = TDv' where td is the average deflection time defined as td = f]mcc{eB)~^ . 

In order to compute a probability Pc{v) for particles having a given velocity v in the downstream 
reference frame, we need to average over all v' velocities corresponding to a a given v downstream. 
In the reference frame of the average flow within the shock transition, velocities v' corresponding 
to the same v are different in different directions, namely v' = v — {{u) — U2). The angle-averaged 
value of v'^ for these particles is given by 

1 —"2 

{{vx - U2{cr - l)/2f + vY) dv^ j/ j dv^ 

— V 

n cr - 1 cr^ - 1 n 

- + ^^vU2 + ^^ul. (9) 

Let us assume that, on the average, electrons in the incoming flow are thermalised to the 
level Xs in -^c collisionless interactions. For the sake of illustration, let us calculate the number 
of scatterings, iVc, which yield a given injection efficiency. The involved time is approximately 
Atx2 = NcTu and therefore 

{u)Ax 3{uf Ati2 _ 3(u)2 
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-^c, (10) 



SO that the probability becomes 



^c(.') = exp(-^i^f^)Vj. (11) 



4 \v' 



Eq. lO, together with probabilities ©, lfTT)l and Eqs.®, shows that, in order to get an 
injection efficiency <r = 10~^, A'c must be equal to 9 for Xs = 0.001, and to 770 for Xs = 0.1. 



It is interesting to note that our expression for Pc behaves Hke the "leakage probabiUty" Vcsc of 
Malkov calculated for protons. Namely, the probability for protons to leak across the shock 

from downstream is approximately lycsdv) oc exp (-const (y') ^) [3. 

Finally, we want to stress that the introduction of the probability Pc does not affect the slope 
of the accelerated spectrum at relativistic energies. Following standard test-particle approach to 
Fermi acceleration a power-law momentum distribution of relativistic particles is generated, 
with an index a = —(2 — <j)l(a— 1) that depends only on the shock compression ratio. This index 
is obtained by combining the term for the momentum increment per cycle (Ap/p) with that for 
the difference (1 — -Pr) per cycle, in the high-velocity limit. The asymptotic behaviour of both 
terms is oc v~^\ while Eq. JZjl is such that (1 — Pc) oc Therefore, in the high velocity limit Pc 
gives negligible contribution to the formation of the particle spectrum in comparison to Pr, and 
does not affect the formula for a. 



3 Thermalisation of electrons and injection 

In the shock-front reference frame, if upstream electrons and ions enter the leading edge of the 
shock transition with the same velocity then the electron energy is lower than that of protons by 
a factor me/wp. Therefore, if the velocities of electrons and ions are randomized independently 
within the shock front, we obtain Tos = (me//imp)Ts ^ Tg (i.e. Xs = Xo ^ 1 in our notation), 
while the ionic temperatures are about (ne -l-ni)Ts/ni (where "i" denotes ions), namely much closer 
to Ts. 

The temperatures of electrons and ions may get closer, if there is a process within the shock 
which allows energy exchanges between the two species. In collisional shocks, the equilibration 
process is Coulomb scattering between ions and electrons, while, in the collisionless case (like it 
generally occurs in SNRs), turbulence plays the dominant role. 



3.1 Results from observations 

Since longtime, it has been suggested that plasma instabilities could lead to prominent heating 
of electrons within the shock (e.g.f2Qj). Some observations and theoretical results put forward 
the possibility that collisionless processes within the shock of SNRs could heat electrons up to 
the level Xs ^ 0.4 ([2] and references therein). Results on SNR DEM L71 in LMC j2il and on 
RCW86 also suggest Xs ~ 0.3. Analysis of Chandra data on Tycho SNR indicates that Xs < 0.1 
dlllHl- Other recent observations (SN1006, Tycho, IE 0102.2-7219) favour a considerably lower 
thermalisation level, namely Xs < 0.03 H- 0.07 [T3 l [2 ^ l9l[T l [T( H [T5l). 

It is important to know how does the level Xs depend on the properties of the shock. Ob- 
servational estimations of the shocks with Mach number up to ~ 400 suggest that stronger 
shocks (namely with higher T4) could equilibrate species less effectively. Namely, Schwartz et 
al. ([21]) present results of measurements of Tca/Tg for interplanetary shocks and planetary bow 
shocks {Ai < 25) and find strong evidence that this ratio depends on the Mach number as Ai^^. 
Ghavamian et al. estimations for a number of SNRs seem to extend this trend to stronger 
shocks, with 25 < M < 200. Rakowski jJH] summarises the observational methods and estima- 
tions of Xs in SNRs shocks and confirms the inverse dependence in the range 25 < M < 400. 



3.2 Results from Bykov and Uvarov (1999) 

Interactions of electrons with ion- or self-generated waves could be responsible for both accelerating 
and heating of electrons (see |19l Hj for a review) . 

BU99 have considered the interactions of electrons with ion-generated electro-magnetic fluctu- 
ations and have developed a kinetic model that accounts at the same time for electron injection, 
acceleration and thermalisation in quasiparallel shocks. Their model is applicable for shocks with 
local Mach number Ai less than ~ y^m^Jrn^. They have introduced the effective electron tem- 
perature Tcff (measured in units of the upstream temperature To), which may be related to our Xs 

by 

To 

^^ = T,^-=Tos—-—-, (12) 
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Figure 1: Injection efficiency <; versus post-shock ratio Xs for a model of BU99 for Bohm-like 
diffusion and diffusion boundary conditions. Curves are labelled by their respective Mach number. 
The approximation = 2 x 10^ [AA'^Xs) is shown by dashed lines. 



and have shown that it depends on the Mach number. This dependence can be approximately 
described by a power law: T^g oc ^A°^, with index < a < 2 depending on which model of 
wave-particle interaction is considered. Therefore, the level of thermalisation depends on the 
velocity of the shock: Xs « M"^^^ for strong shocks, namely the higher the velocity the smaller 
the thermalisation level. 

If A' is momentum independent (transport of electrons is due to large-scale magnetic field 
fluctuations that provide effective heating) then a « 2, and the level of equilibration Xs does 
not depend on the Mach number. This means that Tcs cx; V^, but with a factor that may be 
higher than that inferred from Rankine-Hugoniot equations for the electron population [21]. The 
opposite case is when electron heating in the shock transition region is effectively suppressed by a 
developed small-scale vortex turbulence, giving a ~ 0. In such a situation the postshock electron 
temperature is Tcs ~ To, independently of the Mach number. Another interesting model of wave- 
particle interactions is Bohm-like diffusion, for which X' (x p and a = 1. In the present paper 
we consider only the Bohm-like diffusion case since it seems to be in agreement with observations 
(namely Xs oc , see Sect. l3.1|l . 

BU99 also introduce the dimensionless parameter T = uiAx/v'X' (calculated for electrons 
with V = fth), and in their Fig. 4 they show its dependence on Tcff, for different models of wave- 
particle interactions (Fig. 4a for diffusion boundary conditions and Fig. 4b for free escape boundary 
conditions). In particular, their curve 4 represents results for Bohm-like diffusion. 



3.3 Application of BU99 results to our model 

In our paper we have approximated BU99 numerical results by using r = T^ly^ - 1 where ^ IS a 
constant. This, together with Ijl2ll . gives 

'(« + l)M'x,\"° 1 



In a Bohm-like case, i.e. with a ~ 1, ^ ^ 1.25 corresponds to the diffusive boundary conditions and 
^ = 0.75 to free escape boundary conditions. The parameter T is proportional to the combination 




Figure 2: Electron distribution functions and probabilities calculated for two different values of 
Xs- 1 - Maxwellian distribution, 2 - probability to return to the shock Pr, 3 - probability to 
cross the shock Pc, 4 - initial distribution /inj of electrons injected into acceleration process, 5 
- final hybrid electron distribution {a = 2). Plots are calculated for the model of electron- wave 
interaction developed by BU99 (case of Bohm-like diffusion with diffusive boundary conditions), 
M = 40. a) Xs = 0.01, in this case y^in = 0.12, yb = 2.6, ^ = 0.048, = 20; b) Xs = 0.06, in this 
case 2/min = 0.05, = 3.9, <j = 3.1 x 10"^^ = 930. 



{u)Ax/k in the exponent of the transition probability. This allows us to write 

P( KA\ ( 3(a + i)r(xs,A^) ^ 

P,{v,XsM) = e^V^ ^ ^75 ) (14) 

for nonrelativistic electrons and A' oc p. 

By using Ijl3|l for T in II14|I . the dependence of the fraction of injected particles on the level 
of electron thermalisation Xs may be obtained (Fig.nj. By comparing Ijl4|l with Ijlljl . using ((HJ for 
ymin, we finally obtain a relation between F and Nc'- 

A(a ~ l)xsr(xs) .... 

= 7 — T ■ (15) 

(CT + l)<TXo 

The calculated dependence of the injection efficiency on the Mach number and on the thermali- 
sation level is shown in Fig.n]for the Bohm-like diffusion. The range of values plotted in this figure 
corresponds to a range from 2 to 20 for F in Fig. 4 of BU99, with the minimum F corresponding 
to the maximum The curves in Fig. are essentially the same curve, with different horizontal 
offsets. Namely, the formula for the injection efficiency <;(A^,Xs) is very well approximated by a 
function of Ai^Xs- The reason of this can be found in Eq. (EJ, together with the explicit definitions 
of Pr and Pc (respectively, Eqs. Q and 10): for standard parameter ranges, the most effective 
term is the argument in the exponential of Pc, which is proportional to F. In turn, Eq. (I13II shows 
that the dependence of F from M and Xs is only through the combination M^Xs- In this sense, we 
may say that F is a function of a single parameter (not considering, of course, the dependence on 
the assumed diffusion type and on the boundary conditions, which can be accounted for by using 
parameters a and ^). For instance, in a Bohm-like case a power-law approximation of the curves 
shown in Fig-^is <^ ~ 2 x lO^A^^^^xr^ (this approximation is represented in the figure by dashed 
lines). Since in this case Xs oc M.^^ , the overall dependence of the injection efficiency on the Mach 
number in a Bohm-Hke case is as strong as <x . 

In order to allow for different types of diffusion, as well as for different electron-wave interactions 
etc., one could consider a more general case, in which: i) Xs M~™; and, ii) <r oc M~^''x7'' (where 
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Figure 3: Velocity corresponding to the maximum of the distribution function /inj versus Xs for a 
few choices of the shock Mach number (see curve labels) . Model of electron kinetics is the same 
as in the previous figure. 

we expect q to be always positive). We then obtain cx with b = q{2 — m). In other words, in 
the case of a decelerating SNR shock Xs always increases, while increases if m < 2, and decreases 
if m > 2. 

4 Discussions and Conclusions 

The electron injection and thermaHsation are not independent processes. This is clearly outHned 
by Fig. 121 where the probabilities Pr and Pc as a function of reduced momentum are shown for two 
values of Xs together with the initial distribution I^J of injected particles. The hybrid electron 
distribution ncsfii(y)dy, Maxwellian up to yh and power-law above |22| . is also shown on the 
figure to see the differences. The break momentum yh is given by the assumption that all injected 
particles obtain momenta higher than y^ after acceleration i.e. is defined by = J^^ f^{y)dy. 

It is a common believe that only particles from the energetic tail of MaxweUian distribution 
are capable to be accelerated. On the contrary, the distribution /inj shows that thermal particles 
with velocities v > Umin have the possibility to participate in acceleration process, although with 
different probability. The minimum velocity fmin = 0.07 (xs/0. 03) fth may be considerably less 
than the thermal velocity (see Eq. Q). The most probable velocity at which the maximum of 
the distribution /inj occurs is v^, w (2-;-3)wth for a wide range of injection fractions ^ = lO^'^-^-lO^^ 
(Fig. Eland compare with Fig.QJ. In other words, most of the electrons are injected with velocities 

The injection efficiency of electrons in a collisionless shock is associated to the process of 
electron heating within the shock, through the competition of two effects. On one side, the higher 
the post-shock electron temperature, the higher the energy of thermal electrons and the higher the 
fraction of those which are ready to cross the shock from downstream to upstream (this is given by 
the probability Pr, lines 2 on Fig|2ji,b). On the other side, however, the higher the temperature, 
the higher the number of scattering centers. Electrons traversing the shock from downstream to 
upstream also interact with these sites and the more such interactions the less electrons are able 
to cross the shock and to enter into the Fermi acceleration loop (see probability Pc, lines 3 on 
Fig [2^1, b). In this paper we show that, for a given M, the combined effect of these two processes 
is that the quantity <r decreases with increasing of Xs (Fig.QJ. 

Both injection and thermaHsation are sensitive to the Mach number. It is shown (see Sect. l3.2|l 
that, for a standard range of parameters, <;(xs, A^) is a decreasing function of a single argument 
<j = '^{Ai'^Xs)- Theoretical models show that in high- velocity shocks the energy of the shock is 
transferred to the thermal electrons less efficiently, so that Xs c>c M~"^ with < m < 2 (BU99). 
Observations favour a dependence Xs oc , suggesting a Bohm-Hke type of diffusion. Our 

calculations show that the level of electron-ion equilibration is expected to depend on the injection 
fraction as well, so that the approximate relation between these three parameters is Xs « A4^^<;^^^'' 



(Sect. ESI- The smaller the Mach number, the higher the level of electron-ion equilibration for 
a given injection efficiency. On the other hand, for a given thermalisation level, the stronger the 
shock the less particles can be injected. 

To conclude, we would like to review the assumptions used in the present paper. 

Our approach is in test-particle approximation. Actually, it is known that, in young SNRs, 
shocks could be strongly modified; and the inclusion of nonlinear effects could change our results 
significantly. The usage of a test-particle approach in the present paper is however consistent 
with what done by BU99, and the results of their analysis are valid at least for shocks with 
Alfven Mach number less than ~ 43. Thus our results should be appHcable at least to SNRs 
either in the late adiabatic phase or beyond. Our opinion is that, together with using nonhnear 
treatments, it is valuable investigating what happens in the linear (test-particle) approximation, 
also in consideration that a nonlinear theory has anyway to give, as limit case, the linear results. 

Eqs. (QJ and assume nearly isotropic distribution of particles. In general, this is not fully 
true for the thermal population right after the shock. In order to overcome this difficulty, in 
numerical calculations these formulae are assumed to apply a few mean free paths downstream, in 
order to insure that the distribution is isotropic in the local frame (e.g. jQ). Within our approach, 
this implies some restrictions on the underlying physics. Our assumptions about properties of the 
scattering centers in our model (see Introduction) require that the timescale for isotropisation is 
not larger than the timescale for one interaction. This means that, in order to assume isotropy of 
particles velocities, we would in principle need to increase the number of interactions Nc at least 
by one (see Eq. lfTT|l V Since already iVc >• 1 for, say, Xs > (see estimations after Eq. dJ), 

this increment would not change much our results, in the case of shocks producing an electron 
population thermalised up to the level Xs higher than 10"'^. Since Xs oc Ai~^, our models are 
Hmited again to the shocks with moderate Mach numbers. 
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